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METHOD FOR CHARACTERIZING SAMPLES 




^/^y^ The present invention relates to methods for characterizing samples. These methods are i.a. 
used to investigate test samples from a production, patients or samples collected in any 
5 other way. 
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Background of the invention 

When a sample is to be characterized for components, the components are generally sepa- 
rated from each other in a first step in order to identified and quantified in a later stage, 
1 0 However, it is not always possible to separate the components or it may not be motivated 
from a time/cost benefit reason. The samples may then be characterized spectroscopically 
whereby the components are identified by means of their unique spectral responses. 



If one has a collection of samples and is aware of which components they comprise, it is, as 
15 a rule, trivial to determine their concentrations spectroscopically. This is due even if the 
spectral responses of the components overlaps each other. If, however, the components are 
*p unknown, the problem is muck more complicated . The situation was analysed for the first 

f e time in detail by the mathematics Lawton and Sylvestre ( Technometrics . 13, 617, (1971)), 

fU who showed that it is impossible to find an unique solution even for a 2-component system. 

20 In 1990 we developed an experimental method, which partly solved this problem (Kubi&ta, 
Chemometrics and Intelligent Laboratory Systems. 7, 273, (1990)). We then showed that if 
P one carried out two spectroscopic measurements on each sample, in stead of one as pre- 

viously used, and the measurements were such that the contribution of the components to 
these measurements had the same distribution of the intensities, but of different magnitude, 
25 then both the spectral responses as well as the concentrations of the components could be 
determined Mathematically, these measurements are described using the equations: 

A ™ CV or aj (i)=Xv,W j ~ U ....n 

r 

30 B = CDV or bj (A) = 2 C ^> V / j - 1.2 ....n 



wherein A is a matrix comprising spectra of the first type measured on the n samples; B is a 
matrix comprising spectra of the second type measured on the same n samples; C is a 
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matrix comprising the concentrations of the r different components in the n samples; V is a 
matrix comprising the normalized spectra of the components; and D is a diagonal matrix, 
the r diagonal elements of which being the ratios between the responses of the components 
obtained in the two measurements, All spectra are digitalized in m points. We showed that 
5 the concentrations of the components (C), their normalized spectral responses (V) and the 
ratio between their responses obtained in the two measurements (D) could be determined 
only outgoing from the information obtained from the spectra as measured (A and B). We 
further described how the number of components of the samples (r) could be estimated, 

1 0 One restriction using this method is that the number of components are not allowed to ex- 
ceed the number of samples, which from a practical point of view means that the method 
can not be utilized on smaller series of samples and can not be applied on the whole for 
analysing isolated samples. 

15 Several spectroscopic techniques, such as fluorescence, nmr, etc., can generate 2-dimen- 
sional data described by the equation: 



where the signal, I(at$, is determined as a function of two variables, a and (3, and are the 
20 sum of the contribution of the components in each point, which contribution is proportional 
to their concentrations (Cj) and the products of their (normalized) 1 -dimensional responses, 
l£a) and Out of these responses the components can be identified. In a steady state 
fluorescence spectroscopy//^ and I/fi) are the excitation- and emissions spectra of the 
components and are, as a rule, designated I?(AJ and I ( cm (A er J, wherein X tx and X cm are the 
25 excitation and emission wavelengths. The shape of an excitation spectra of a pure com- 
pound is, in general independent of the emission wavelength used at the measurement, and 
the corresponding is due for its emission spectrum. The fluorescence signal monitored, if 
necessary after a correction for the inner filter effect (Kubista et al, The Analyst. 119, 417 
(1994)) ? is proportional to the concentration of the compound. In a sample containing more 
30 compounds the total signal is the sum of the contribution by each component. As fluore- 
sccnce is measured in an arbitrary unit, eq. 1 contains a proportionality constant (k)« 



r 



The information of the 2-dimensional spectrum l(ccp) is insufficient to unambiguously 
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determine the spectral responses of the components. Different approximative ways have 
been suggested but these do not function sufficient satisfactorily even for a 2-component 
mixture (Burdick and Tu, J. Chqnometrics. 3, 431, (1989)). 

5 The present invention is a method for analysing isolated test samples, or a couple of test 
samples without using references in such a way that the components can be identified. 

description of the figures. 

Figure 1 . Emission spectra monitored using different excitation wavelengths using a parallel 
10 polarized light (above, left) and a perpendicularly polarized light (above right), respectively, 
Down to the left the calculated emission spectra of the components are shown, and down to 
the right the calculated excitation spectra of the components are shown. 



\i Figure 2- A) Excitation spectra registered using different emission wavelengths from two 

Ui 15 solutions containing POPOP, dimethyl POPOP, antracene, and diphenyl antracene* B) The 

Ql excitation spectra of the components as calculated. 

Figure 3. A) Emission spectra registered using different excitation wavelengths of two 
p solutions containing POPOP, dimethyl POPOP, antracene, and diphenyl antracene. B) The 

[r% 20 excitation spectra of the components as calculated. 



D 



Brief description of the invention 

The present invention is a method for analyzing test samples in such a way that its compo- 
nents can be identified without the need for any reference data. The method is based upon 
25 the following four steps: 

1, The test sample is analyzed using a method generating a 3-dimensional response accord- 

ingto: itoM-ti^efiWir), 

wherein r is the number of components contributing to the signal, and 1/ cfy and Irffi) and 
30 %{y) are the arbitrarily normalized I -dimensional responses of the components, which 
responses normally consist of spectral or concentration variations. 

2, The number of components r as the samples contain is estimated, 

3, For each component its 1 -dimensional responses l t ($ and 10 and // y) are determined. 
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4. Out of the responses, the components are identified. 
Defojied description of the present invention 

As the title indicates the present invention relates to a method for characterizing isolated test 
samples in a way that makes it possible to identify its components without any need for 
using reference samples. This is done through a strategic design of experiments which 
makes it possible to register a 3-dimensional response being proportional to the concentra- 
tions of the components, and the contribution from each component is the product of its 
specific 1 -dimensional responses: 

/(flr,Ar)-=S^(«y,W,(r) 

Such registration can be carried using certain forms of fluorescence spectroscopy, e.g., by 
means of a time disintegrated monitoring of emission/excitation spectra, i.e., the signal is 



registered as a function of excitation wavelength, emission wavelength, and time 

y 

f-i 



15 /(^,^,rt-£c,/,<^)/,(^)J,(0 

In these cases it is often suitable to gather the concentration of the components ci and the 
time declinations to a time dependent concentration: 



FU 

*d 20 The time can be time after light pulse (whereby Cj(t) is proportional to the fluorescence 

declination), time after mixing of e.g., a stop-flow experiment (whereby Cj(t) is the variation 
of the concentration of component i with time), time after treatment, such a photo bleaching 
(selective destruction of certain components using light), chromatographic or other form of 
separation, etc. At the analysis of such data the concentration variation of the components 
25 are calculated, as well as their excitation and emission spectra. It is of interest to note that 
intermediate components which are neither present at the beginning (Cj(0) = 0) or at the end 
( c i(") = 0) of the experiment can be identified from its calculated spectra. 



30 



There is a further possibility in varying the polarization of the light: 
1 (^>^ cmi ^) = (X ex )I i {A m )I t {a) 

or, if the phase-modulated light is utilized, the frequency of the modulation: 
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etc. 

There is further a possibility in varying the outer parameters which influences the concentra- 
tions of the components, such as temperature (pressure, volume, etc.): 

or outer parameters which influence the intensity of the responses of the components, such 
as external magnetic fields (electrical fields, etc.): 

K^„,M) - J>,/, U m )I, (M) 



10 



The spectroscopic technique need not be a fluorescence technique. The method can be car- 
ried out using most techniques which generates 3-dimensional responses, e.g., nuclear 
magnetic resonance spectrometry (NMR) mass spectrometry, etc. It can further be carried 
out using most techniques generating 2-dimensional responses if the responses of the 
1 5 components influence external parameters. Finally, the method can be used using a tech- 
nique generating 1 -dimensional responses, as well, but then it is necessary that two external 
parameters are varied simultaneously and that their influence on the responses of the com- 
ponents are independent so that their contribution can be factorized, 

20 The invention requests that at least two data points are determined in each of the 3 dimen- 
sions, i.e.: 

l^a) wherein a n 1>2 

//^wherein p„(J 2 p m m>2 

I/y) wherein y 1^2 Yia>2 

25 

To determine two data points only in all dimensions are, however, of particular meaning as 
the tolerance of the responses calculated then as a rule is insufficient to be able to identify 
the components. On the contrary it is quite excellent to have two data points only in one of 
the dimensions , e.g, 1=2 (and m»2, and n»2). This exhibits the advantage that the nume- 
30 rical treatment of data is made easier as the responses of the components can be calculated 
using fast algorithms such as Procrustes rotation and GRAM (Kubista, Oiemomefrics and 
Intelligent Laboratory Systems. 7, 273, (1990); Wilson, Sanches & Kowalski, J. Osmo- 
metries, 3, 493, (1989)). In the general case when all 1, n, and m are greater than 2, the 
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solution method is much more complicated and thus considerably more time consuming 
(Liwo, et al, Com puters Chenu 21, 89-91, (1997)). Furthermore, it is quite often of interest 
to carry out the experiment in such a way that one of m and n are considerably greater than 
the other. The reason hereto is that it as a rule, is sufficient for the identification of the com- 
5 ponents, to determine one of their 1 -dimensional responses with a high accuracy. 

The invention is not limited to determinations that generates 3-dimensional responses but 
even responses of a higher order can be used. In general it should be satisfying that the 
response is linear and that the contribution from each component shall be the product of its 
1 0 1 -dimensional responses : 

□ /(*, frr.s....) = 5>,/. (ay i <#/,. (/)',- <*)...- 

m 

^ Of course, the higher the dimension is the more time consuming the numerical treatment of 

1 5 the determined data will become. However, with regard to the very fast development within 

B[j the computer area this will hardly be a practical limitation in the future. 

Si 

The samples to be analysed shall contain substantially the same components, and these shall 

ft! , . . . n 

'i2 be present in different, relative concentrations. The samples are analysed m pair using a 2- 

t-Q 20 dimensional method which provides a response which is proportional to the concentrations 

of the components and the product of the 1 -dimensional responses. This can be expressed 

as* 

25 m 

wherein ?(&) and are spectra of the two samples which in the following will be called 
A and B, determined as a function of the variables a and p, r is the total number of compo- 
nents contributing to the spectra, I {a) and are the normalized 1 -dimensional responses 
of the components, and c, A and are their concentrations) respectively. In a steady-state 
30 fluorescence spectroscopy are the normalized excitation spectra of the components, 
I^fAj, and 10 are their normalized emission spectra, ir(^J- 



ffi 
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The information in these spectra is treated in two steps. First the number of components, r, 
5 is determined, and then the I -dimensional responses of the components* 
/*W.andiT^ 

When r has been determined, the spectral responses of the components. 
The equations: 



10 



can be written in matrix form as: 
SJ A=XC A M 
W 15 B=XC»M 

CP wherein A and B are matrixes comprising the spectra determined, X is a matrix comprising 

the nonnalized excitation spectra of the components, M is a matrix comprising their norma- 

s 

lized emission spectra, and C* and C B are diagonal matrixes comprising the concentraUons 
PJ of the components. By renormalizing one of X or M, the equation system can be rewritten 

,f* 20 as: 

O A=XM 

Q 



B-XDM 

wherein D is a diagonal matrix comprising the ratios between the concentrations of the 
components (D*=C*/C A ). Using A and B, X, M and D can be calculated using known 
25 methods such as Procrustes rotation (Kubista, Chemometri cs and Intelligent Laboratory 
Systems. 7, 273, (1990); and GRAM (Wilson, Sanches & Kowalski, J. ChernoTpetrics, 3 S 
493,(1989). 

As a summary, the present invention relates to a method for experimentally studying two 
30 samples spectroscopically so that the information present in the experimental spectra is 
sufficient to determine the number of components of the samples (r), their spectral re- 
sponses of the 1* dimension, their spectral responses of the 2 nd dimension, and 
the ratios between their concentrations (c, A /c*). 
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The most apparent use of the invention is for the analysis of two samples containing 
common components. All components need not be common, but the majority of those con- 
tributing spectroscopically should be in common (Booksh & Kowalski, J. Chemotrics, 8, 
287, (1994)). The number of components is arbitrary and can exceed 2. 

5 

Another use of the invention is to characterize single samples by first dividing them into 
two part samples containing the ingoing components in different proportions. This can be 
accomplished in several ways, e.g., by filtering, extracting, chromatographying dialysing, 
centrifuging, precipitating, splitting the sample by means of an electrical field, etc. Alter- 
10 natively, the original sample can be used as one sample, and an aliquot thereof, which is 
created in such a way that the components are present in other proportions, is used as the 
second sample. This aliquot can be obtained by selectively eliminating certain components, 
e,g*> by means of adsorption, precipitation, freezing, distillation, selective decomposing 
Nf (e. g , y by light, beat, radio lysis), etc. Another possibility is to create two samples from one, 

yjj 1 5 is to change the conditions for the determination, e.g., by changing the temperature, pres- 

ffj sure, etc. Separation methods, such as different types of chromatography are of interest, as 

$ the components are separated in space, and one, principally arbitrary number of samples can 

be obtained which can be analysed in pair. Using spectroscopic techniques which generates 
y, 2 -dimensional spectra in a fast way, then, furthermore, the detection can be made on-line. 

*0 20 

p 

j*3 Another use of the invention is to determine the concentrations of the components m one 

test sample in relation to a standard sample with a high degree of accuracy. The standard 
sample and the test sample are analysed as a pair, and the ratio between the concentrations 
of the components is obtained as the diagonal element of the D matrix. 

25 

2-dimensional spectra wherein one of the dimensions is time, are of particular interest, 
whereby time is related to time after a disturbance such as a relaxation time. Today, there 
are e.g., fluorescence instruments by means of which one can determine complete spectra as 
a function of time after lightening (either directly after lightening using a light pulse, or 
30 indirectly using phase modulation technique). This gives using a as time, and p as wave 
length, the equation system: 



u 

if* 
hi 
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from which r 9 1ft), 1/A) »d (c^/e,*) can be determined. 
5 Example 

The invention will be further illustrated in four examples. 
Example 1 

A sample is characterized using fluorescence spectroscopy, where excitation wave length, 
10 emission wave length, and light polarization are varied (Figure 1), This gives raise to a 3- 
dimensional spectrum according to: 

lW«a m >a) = 1>,/i <4J'i (*) 



In the example 650 different emission wave lengths (m), 1 1 different excitation wave 
1 5 lengths (n) and 2 different polarizations (a = Oo, called parallel polarization, and a - 90o f 
called perpendicular polarization) (1), arc used. From the response determined, In(A„, ^ 
^ c$ , first the number of components (r) is estimated to 2 (using a statistic test and a visual 

inspection of the principal components). Then the component specific responses are cal- 
culated. For this purpose one uses the fact that only two data points were registered in one 
[g 20 of the dimensions (polarization) and rewrote the 3-dimensional response to two 2-dimen- 

0 sional responses. 

1(^,^=0°) = (A m y, W m )J,. (a = 0°) 



m 



- 90°) = |> ( ./,. {A„ )/,• )/, (a - 90°) 
25 These can be described using the equation system: 
I°=Xa°M 
r=Xa M M 

which can be solved using Procrustes rotation (Kubista, Chemometrics and Intelligent labo- 
ratory Systems. 7, 273, (1990)). This gave the normalized excitation intensities of the com- 
30 ponents as matrix XU/AJ), (shown down to the right in Figure 1 ), the normalized emission 
intensities of the components as matrix M (I /(A J) (shown down to the left in Figure 1), and 
the ratios between the responses of components to light of different polarization 
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From the calculated component specific responses, in particular the emission spectra, the 
component could be identified as p-bis[2-(5-phenyloxazolyl)]-ben2ene (POPOP), and antra- 
cene. Finally, by comparing standard spectra of POPOP and antracene the concentrations 
could be estimated to some micro molars. 

S 

Two solutions containing the dye compounds POPOP, dimethyl POPOP, antracene and 
diphenyl antracene in different proportions were prepared. On these fluorescence excitation 
spectra were monitored at several emission wave lengths. The number of components were 
10 determined to 4 using a statistic test, and the excitation spectra of the components (Figure 
1), emission intensities and the relation between their concentrations in the two samples 
were calculated. 

Example 3 

15 On the same solutions as in Example 1 the fluorescence emission spectra were monitored 
using a number of excitation wave lengths . The number of components was determined to 4 
using a statistic test, and the emission spectra of the components (Figure 2), excitation in- 
tensities and the relation between their concentrations in the two samples were determined. 



^ 20 Example 4 

CO 

£3 



Characterization of samples containing the dye compound thiazole orange and the polymer 
poly(dG) was made. The samples were analysed in pairs using 2-dimensional fluorescence 
spectroscopy. They contains thiazole orange and poiy(dG) in the relation 
[thiazole orange]/[poly(dG)] of 0.05 and 0.025. Neither poly(dG) nor the dye compound is 
25 fluorescentic as such but the fluorescence arises when thiazole orange binds to the polymer. 
The samples were analysed in two different ways. In one analysis, the fluorescence exci- 
tation spectra were monitored at different emission wave lengths. The number of fluorescent 
components were identified to two using statistic tests, and their excitation spectra and 
emission intensities were calculated.. In the second analysis, the fluorescence emission 
30 spectra were monitored using a number of excitation wave lengths. Once again the number 
of components was identified to two, and their emission spectra and excitation intensities 
were calculated. 



